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We present new measurements of the scintillation and ionization yields in liquid xenon for low 
energy electronic (about 3-7keV ee ) and nuclear recoils (about 8-20keV nr ) at different drift fields 
from 236V/cm to 3.93kV/cm, using a three-dimensional sensitive liquid xenon time projection 
chamber with high energy and position resolutions. Our measurement of signal responses to 
nuclear recoils agrees with predictions from the NEST model. However, our measured ionization 
(scintillation) yields for electronic recoils are consistently higher (lower) than those from the NEST 
model by about 5e _ /keV ee (ph/keV e e) at all scanned drift fields. New recombination parameters 
based on the Thomas-Imel box model are derived from our data. Given the lack of precise 
measurement of scintillation and ionization yields for low energy electronic recoils in liquid xenon 
previously, our new measurement provides so far the best available data covering low energy region 
at different drift fields for liquid xenon detectors relevant to dark matter searches. 

PACS numbers: 


I. INTRODUCTION 

Particle detection technology based on liquid xenon 
(LXe) has been developed extensively in the last decade, 
thanks to the advance and promise of the direct dark 
matter experiments HHZ!- Among these, the dual 
phase xenon technique has made significant progress 
due to its capability for background identification and 
suppression, and scalability to a large ton-scale target 
mass. Such a technique allows the detection of low energy 
electronic and nuclear recoils down to sub-keV with 
both scintillation and ionization signals. Understanding 
the response of LXe to low energy events becomes 
increasingly important in order to properly assess the 
background responses in the detector and to precisely 
extract the dark matter parameters with a positive 
detection. 

In the last few years, low energy nuclear recoils (NRs) 
were measured extensively either by tagging elastically 
scattered neutrons from a fixed energy neutron source, 
e.g. a DD generator mm, or by modelling the response 
and comparing with a neutron source, usually 252 Cf 
and ArnBe, with a spread of neutron energy at the 
MeV level |TOll4| . Such measurements have yielded 
better understanding of the scintillation and ionization 
properties in LXe for NRs below 10 keV nr (keV nT . denotes 
the nuclear recoil energy, while keV ee denotes the 
electron equivalent energy.). The uncertainty of WIMP 
detection sensitivity was being reduced continuously, 
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especially for low mass WIMPs below 10 GeV/c 2 , due 
to a better understanding of such quantities. 

On the other hand, the scintillation and ionization of 
low energy electron recoils (ERs) were not well measured, 
partly due to the difficulty to introduce such a low 
energy ERs in the LXe target. The Columbia and 
Zurich groups n~5i nr] measured the Compton scattered 
electrons by tagging the scattered gammas at a certain 
angle. However, these measurements are limited to 
scintillation yield with large uncertainties. Recently, an 
attempt to measure ionization yield at 2.82keV ee has 
been realised by using 37 Ar source doped in LXe HU¬ 
SO far, that is the only measurement of ionization 
yield below 10keV ee and it’s only at one given field 
(3.75kV/cm), and also with large uncertainty. 

Understanding LXe’s response to low energy ERs 
are not only important for fully understanding the 
background for which the ERs are the dominant 
contribution so far [dj, but also relevant to extracting 
information from dark matter candidates that produce 
ERs, such as from the axioelectric effect [T8]. Several 
dark matter detectors based on LXe are operated at 
different electric fields, with different configurations 
of light and charge detection. For example, LUX 
detector [6j is operated at a relatively low drift field 
(180V/cm), while the ZEPLIN-III Q3] detector was 
operated at a high drift field (3.9kV/cm). Thus a precise 
measurement of the scintillation and ionization in LXe 
for low energy ERs at different fields become extremely 
demanding. Also the understanding of the low energy 
ERs can provide a basis for the background modeling of 
the solar neutrino in the next-generation large-scale LXe 
dark matter detectors mm- 

Here we report the measurement of ionization and 
scintillation yield using Compton scattered low energy 
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electrons in LXe at different fields. Unlike the 
measurement performed at m fl6] . we did not use 
tagged Compton gammas for a fixed energy ER. Instead, 
we took data for both ionization and scintillation for 
Compton electrons at all different energies, and extract 
the ionization and scintillation yields based on the 
Thomas-Imel recombination model [21] (energy below 
about 7keV ee ). By comparing to the simulation of 
the signal response, we found our measured results are 
quite different from those in the NEST model [221424] 
for ERs, while the NRs give consistent results. We 
also provide the model-independent photon yields for low 
energy ERs with energy from about 3 to 20keV ee , based 
on a combined energy scale. 
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II. EXPERIMENTAL SETUP 

The measurement was performed in a two-phase LXe 
time projection chamber (TPC), with four Hamamatsu 
R8520 photomultiplier tubes (PMTs) on the top and one 
R11410 PMT on the bottom viewing a 1-cm thick LXe 
target, allowing simultaneous measurement of both the 
scintillation (SI) and ionization (S2) signals. The X-Y 
positions of events are reconstructed through the S2 hit 
patterns on the four top PMTs. More details of the 
setup was described in [25]. Low energy ERs in LXe are 
obtained using an external 13 7 Cs source. 662 keV gamma 
rays from the source also provide energy calibration and 
stability monitoring. A 252 Cf neutron source is used to 
produce low energy NRs in the detector. 

During the operation, the anode electrode was 
connected to ground. The gate grid was fixed at -4kV, 
providing sufficient extraction and gas amplification 
fields of about llkV/cm (with the liquid level of about 
2.9 mm) for the electron emission into the gas phase. 
The cathode was adjusted accordingly to provide nominal 
drift fields from 200V/cm to 2kV/cm (correspond to 
236V/cm to 1.92kV/cm according to [25]) across the 
1-cm drift region. Due to the limitation of electron 
transmission through the gate grid, we conducted a 
special high drift field run by lowering the liquid surface 
below the gate grid. During such a run, the gate 
grid was connected to ground and the cathode was 
set at -5kV. Using the mean drift time (Fig. [l]) of 
photoelectrons emitted from the cathode in the form of 
S2 after-pulses [271429] and the electron drift velocities in 
LXe [3Tj, we calculate the total drift length and the drift 
field in LXe to be 8.1±0.5mm and 3.93±0.15 kV/cm, 
respectively, based on the field simulation using the finite 
element analysis software package COMSOL [25]. 

The low energy events relevant to our study are below 
7keV ee , while the energy calibration is based on the 
662keV ee total absorption events. The low energy data 
and the calibration data were taken with different PMT 
gains to avoid saturation for the calibration events. The 
calibration data were taken with the gains of four top 
R8520 PMTs at 2xl0 6 and the bottom R11410 PMT 


FIG. 1: Distribution of time differences between a main S2 
and an after-pulse S2 in the special high field run. The peak 
corresponds to the drift time of photoelectrons emitted by the 
cathode. The mean drift time from the cathode to the liquid 
surface is 3.37±0.19ps. 

at 2.5xl0 4 . For the low energy data, the gain of the 
bottom R11410 PMT was at 4x 10 6 . All SI and S2 signals 
are calculated in unit of photoelectron (PE) by dividing 
the signal outputs by the PMT gains. The signal from 
bottom PMT was split into two channels, one for SI with 
an external x8 amplifier (CAEN N979B) and another for 
S2 directly, and fed into flash analog-to-digital converter 
(FADC) digitizers. S2 from the bottom PMT only 
was used in the analysis due to the high collection 
efficiency and resolution of the R11410 PMT. S2 from 
the bottom PMT was used for trigger and a trigger 
threshold of 58 PE, corresponding to about 1.8 extracted 
electrons, was obtained. In addition, a high-energy veto 
with the upperlimit of about 20000 PE on the S2 signal 
was implemented in the trigger of the low-energy recoil 
measurements in order to assure a good linearity of the 
S2 signal. 

The gains of PMTs used in our study were measured 
every week to ensure PMT stability. The PMT gain 
was first measured at a high reference voltage using a 
weak LED light generating single photoelectrons. Then 
the LED light was increased and signal dependence on 
the PMT voltage is measured, covering the entire range 
of gains used in our measurement. During the gain 
calibration, the rate of LED pulses was set to about 
100 Hz. The gain dependence on the PMT voltage can 
be described by the following equation. 

G(V)=G(V ref )(^j , ( 1 ) 

where V and G(V) represent the PMT voltage and the 
corresponding gain. V re f is the reference voltage, which 
is 800 V for R8520 PMTs and 1500 V for R11410 PMT. 
For the bottom R11410 PMT, the combined calibration 
data during the six weeks of operation show an averaged 
G(V re f) of 3.46x10® and n of 8.33, with standard 




3 


deviations of O.llxlO 6 and 0.09, respectively. The gain 
variation due to different event rate in our measurement 
shall be small (<5%) [303■ Also the quantum efficiency 
(QE) of R11410 PMT shall not depend on the bias 
voltage according to Hamamatsu corporation. 

III. DETECTOR CALIBRATION 

The measurements of the low energy ERs and NRs 
were taken within about two months (Nov. 2013-Jan. 
2014), during which daily gamma calibrations using the 
137 Cs source at a drift field of 987V/cm were carried 
out to monitor the stability of SI and S2 signals. The 
event rate in the daily monitoring data is about 400 Hz 
(with about 30Hz of the triggers due to the background). 
Events in the central region with reconstructed radius 
less than 10 mm and drift time between 3.5-8.5 mm were 
selected to reject regions with bad field uniformity. We fit 
the SI and S2 spectra by an exponential function plus a 
Gaussian, to obtain the SI and S2 yields for the 662 keV 
gamma rays. 



FIG. 2: The evolution of the SI light yield (upper) and 
S2 yield (lower) for 662 keV gamma rays under a drift field 
of 987V/cm. Typical SI and S2 spectra, with fits using 
an exponential function plus a Gaussian, are shown in the 
insets respectively. The decreasing of the S2 yield over time 
is modelled by a linear function with a decreasing rate of 
3.1±0.2 PE/keV ee /day with the S2 yield at the beginning of 
the measurement ( 15 th Nov.) being 1588PE/keV ee . 

Fig. m shows the time stability of SI and S2 yields 
during the two months of operation. The SI light yield 
is stable during the whole period with an average value 
of 3.94PE/keV ee and a standard deviation (s.d.) of 
0.06 PE/keV e e. We observe a decrease of S2 yield over 
time, which is fitted by a linear function with an average 
value of 1514±22 PE/keV ee and a decreasing constant of 
3.1±0.2 PE/keVee/day. The time dependence of the S2 
yields is caused by a micro leak in the system gradually 
lowering the liquid level, reducing the gas field. S2 signals 
of the low energy measurements are corrected for such 


a time dependence according to the linear fit shown in 
Fig. [2] 

The detected SI and S2 signals can be written as SI 
= PDE- N 7 and S2 = EAF- N e , where PDE and EAF 
represent the photon detection efficiency and electron 
amplification factor, respectively. PDE is the product 
of the light collection efficiency, which is related to 
the detector geometry, and the quantum efficiency of 
the PMTs. EAF is the product of electron extraction 
efficiency and the gas gain. N 7 and N e are the number 
of scintillation photons and drifting electrons after the 
electron-ion recombination process, thus PDE and EAF 
are independent of drift fields. The PDE is mainly 
relevant to the liquid level because of the total reflection 
of the scintillation light on the liquid-gas surface. The 
EAF mostly depends on the gas held strength and the 
thickness of the gas gap. The photon and electron yield 
at various drift fields were measured extensively before 
and their values at 987V/cm are 25.43±1.02 ph/keV ee 
and 47.56±1.90e _ /keVee, respectively, from NEST 
v0.98 [52]. Based on these values, we obtain a PDE 
of 15.5±0.2% and an EAF of 31.8T0.5 PE/e" from 
our calibration data. The systematic uncertainties for 
the PDE and EAF are estimated to be ±1.3% and 
±2.6PE/e _ , respectively, which take into account the 
global uncertainty of 4% for the NEST prediction, 
the S2 yield uncertainty of 1.2% induced by 20% 
electron lifetime variation during the operation, and 
the gain difference uncertainty of 7.2% between the 
monitoring 137 Cs calibration and the low-energy recoil 
measurements. 

For the special run with drift held at 3.93kV/cm, 
because the liquid surface was adjusted below the gate 
grid, the PDE and EAF are different from the normal 
runs. The 137 Cs calibration during the special run gives 
an average SI light yield of 3.99±0.03 PE/keV ee and an 
S2 yield of 1053±5 PE/keV ee . The photon yield and 
electron yield at 3.93±0.15 kV/cm for 662 keV gamma 
rays are 19.99±0.81 ph/keV ee and 52.67±2.11 e - /keV ee , 
respectively [52]. This leads to a PDE of 20.0±1.7% 
and an EAF of 20.0±1.7PE/e~. The S2 signal’s time 
dependence is negligible because it took only two days 
for the special run and the detector is stable within this 
time scale. 


IV. RESULTS AND DISCUSSION 

Low energy ER data were taken at different drift helds 
from 236V/cm to 3.93 kV/cm with a low-energy event 
rate of about 30 Hz. To accumulate enough statistics, 
each measurement was taken for about 24 hours. We 
also took the NR data at these drift helds with an event 
rate of about 10 Hz. Single scatter events in the same 
fiducial volume as for the calibration data are selected 
for the analysis. As an example, the low energy ER and 
NR bands at a drift held of 236 V/cm are shown in Fig. [3] 
In order to extract the photon and electron yields from 
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FIG. 3: The measured (left) low energy ER and NR bands under a field of 236V/cm, along with the bands from simulation 
(right) taking into account the detection efficiencies (PDE and EAF) and statistical smearing effect. The magenta (blue) and 
violet (red) solid (dashed) lines are the means of the ER and NR bands, respectively, in data (MC). The input of scintillation 
and ionization yields to the simulation for the NRs is based on the NEST VI.0 [231 [24 :. The input for the ERs is based on a 
y 2 analysis as discussed in the text. 



FIG. 4: The ER (blue) and NR (red) spectra obtained 
through the Geant4 simulation. The spectra with energy 
below 20 keV e e and 50 keV nr are used in the signal band 
simulations for the ER and NR, respectively. 

our measured ER and NR bands, we consider the signal 
generation mechanism in LXe [33]. 


e = (iV 7 + Ne)Wq 

(2) 

% = {N ex /N i +r)N i 

(3) 

N e = (1 - r)Ni 

(4) 


where e is the energy deposition, W q is the average 
energy required to produce a quanta (photon or electron). 
W q is found to be 13.7±0.2eV [22], independent of 
energy deposition and drift field. N ex /Ni is the ratio 
of number of excimers formed N ex to number of ions iV) 
created from the energy deposition. N ex /Ni is taken as 
a constant at about 0.06 ;22j for ERs and is modelled as 
a function of the applied field and deposited energy for 



FIG. 5: The best-fit values of 4 £/Ni obtained by comparing 
Monte Carlo (MC) and data for ER. The inset shows the y 2 of 
the scanned 4£/lVi to the ER means in 236 V/cm data, which 
obtains a minimum x 2 = 16 (number of degrees of freedom = 
7) at 4£/lVi=0.0214 with a statistical uncertainty of 0.0003 
(with Ay 2 —1). The red dashed line represents the prediction 
by NEST. The black rectangles are from E. Dahl [3Ji (ER 
only data). The blue shadows represent the systematic 
uncertainties, which are induced by the uncertainties of the 
PDE and EAF. A fit through our data points gives 4 £/Ni = 
0.046^o ooe] E -0 ' 140 . The index of -0.14 is noticeably close to 
the calculated value of -0.1 by Eric Dahl [34] and the NEST 
value of -0.11 [22] . 

NRs [23]. r is the electron-ion recombination fraction. In 
the Thomas-Imel box model [21] approximation for low 
energy events (which is used in NEST [22]). 

(5) 

Here /it and E are the mobility in xenon and the field 
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ER means (data) 

NR means (data) 

ER means (MC using NEST scale) 

NR means (MC using NEST scale) 

ER means (MC using scale in this work) 


FIG. 6: The measured mean values of ER (blue data points) and NR (red data points) from Logio(S2/SI) bands for drift fields 
from 236 V/cm to 3.93 kV/cm. The red dashed lines are the NR band mean lines from simulation using NEST model |231I24| . 
The blue dashed lines are the ER band mean lines from simulation using the photon and electron response obtained in this 
work, as shown in Fig. [5] For comparison, the ER band mean lines from simulation using the NEST model are shown as the 
green dashed lines, which are consistently lower than the measured values from this work. 


strength, respectively, a and a are the recombination 
coefficient and the box volume size ED- The parameter 
4£/iVj = a/(a 2 fj,E) is a dimensionless constant at a given 
drift field [251 EH] • 


For a given recombination fraction r, the photon and 
electron yields at a given energy can be obtained based 
on Eqs. SH By taking into account the detector related 
parameters, PDE and EAF, and the statistical effects, we 







































6 


simulate the ER and NR bands for comparison with the 
data bands, as shown in Fig. [3] (right). The simulation 
takes into account the Poisson fluctuation of photon 
detection and Binomial fluctuation of the electron-ion 
recombination, as well as a Gaussian fluctuation on 
the recombination fraction Ar(e). Additionally, the 
fluctuations caused by the PMT’s single-photoelectron 
(SPE) resolution and the gas gain, which is not significant 
due to a relatively large number of electrons, are taken 
into account. Both the electron and NR energy spectra 
are obtained from a Geant4 simulation jJ5]. For the ER 
band, the simulated spectrum below 20keV ee is used. 
The ER and NR spectra from simulation are shown in 
Fig. [4] 

A y 2 analysis to compare the measured and simulated 
ER band means is carried out by scanning different 
4£/Ni. The minimized % 2 value corresponds to the 
best-fit 4 £/Ni value (see Fig. [5] inset). The SI range 
is from 8 to 40 PE. The 8 PE threshold is constrained 
by the high dark rate observed in data. The SI 
range corresponds to the energy range of about 3 to 
7keV ee , based on the best-fit 4£/Ni in this work. 
The recombination fluctuation A r can be obtained by 
comparing the band widths between data and simulation. 
It is out of the scope of this paper and will be reported 
later. The uncertainty of Ar brings an uncertainty 
of io 002 the MC ER means, which will be taken 
into account in the interpretation of the uncertainties of 
4£/N 4 . 

The most probable values for 4£/A/ obtained from 
the x' 2 analysis for all scanned fields are shown in 
Fig. § The systematic uncertainties are dominated by 
the uncertainties of PDE and EAF. A fit through the 
ER means under 236 V/cm using Eq. [5] gives a value of 
0.0214±0.0003 for 4£/lVj, which is ~30% lower than the 
value in NEST. 

Our measured 4^/ATj values are significantly lower than 
those measured by Dahl |34|. The derived predictions 
from NEST, which uses mainly the data from Dahl, thus 
give a higher value than our measurement. It is because 
that Dahl’s measurement was from a LXe detector 
without X-Y position sensitivity and thus the edge effect 
gave a large systematic error. In our measurements, 
the ER photon yields from data in the entire volume 
without radius selection are observed to be closer to 
the NEST predictions, although still about 2.5ph/keV ee 
lower than the NEST predictions. The ER photon yields 
derived from data in the central fiducial volume have 
larger difference, which is about 5ph/keV ee lower than 
the NEST predictions. 

Fig. [6] shows the measured ER and NR band means 
at all scanned fields, together with the best fit curves 
from simulation. For data at each field, we also plot 
the simulated band means using the parameters in the 
NEST model [23]. The ER band means from the 
NEST prediction are consistently lower than from our 
measurement. For NRs, the band means from NEST 
and our measurement agree very well. 



FIG. 7: The comparison of photon and electron yields 
derived from our measurement and NEST to the fixed 
energy measurements in HZI and }36j . The red circle is the 
measured ionization yield of 2.82 keV ee gamma under a field 
of 3.75kV/cm |T7], and the blue rectangles are the measured 
photon yields under a field of 450V/cm [U51 . The blue 
shadows represent the systematic uncertainties of the photon 
yields from m■ The statistical and systematic uncertainties 
of the photon yield from j3Sl are the propagations of the 
uncertainties of the relative yield R e and the SI quenching 
under 450V/cm q(450). The red (violet) shadow represents 
the range of the charge (photon) yields at 3.75kV/cm 
(450V/cm) based on the Thomas-Imel box model and the 
updated 4£/Ni obtained in this work with the uncertainties 

(Fig- 

Using the measured 4£/ATj and its dependence on the 
field, we are able to predict the photon and electron yields 
at any given field from 236V/cm to 3.93kV/cm for low 
energy electron recoils below 7keV ee - There were very 
few measurements of photon and electron yields with a 
fixed energy source. The only measured electron yield is 
for 2.82keVee ERs at 3.75kV/cm [T7]. For photon yields 
below 10 keV ee with a drift field, the only measurement 
is from the tagged Compton scattering experiment at 
450V/cm [3H]. Fig. [7] shows the predictions using 4£/IVj 
values from our measurement and NEST, compared with 
the fixed energy measurements. For the photon yields at 
1.5, 2.6, 7.8keV ee events in [3B|, we derive the values as 
the product of their relative light yield R e to 32.1 keV ee 
events from 83m Kr at zero field, their SI quenching 
q(450) under 450V/cm and NEST predicted photon 
yield at 32.1 keV ee , at which energy much more accurate 
measurements are available, providing more precise 
predictions from NEST. Due to the large uncertainties 
associated with the fixed energy measurement, the 
ER yields from our measurement and NEST are both 
compatible within the errors. Further measurement from 
other groups using the band comparison method used 
here or more precise measurement with fixed energy 
sources are needed to reconcile the differences between 
this result and those from NEST and Dahl. 
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FIG. 8: (Top panel) The best-fit values of N ea ,/Ni obtained by 
comparing ER means in MC and data. The 4£/Ni and N e2! /N, 
are both set to be free parameters in the minimum-^ 2 fitting. 
The mean N ea ,/Ni is 0.11 with variance of 0.07 among different 
fields. With the N e:c /Ni=0.11, the 4£/Ni obtained in this 
work changes to 4£/Ni=0.042|/°'°°gj E -0 ' 140 correspondingly. 
(Bottom panel) The simulated ER band with N ea ,/Ni=0.11 
and 4£/Ni = 0.0197 (compatible with the ER measured under 
236 V/cm in this work). The red (black) solid line represents 
the mean of simulated ER band with N e;E /Ni and 4£/N; of 
0.06 (0.11) and 0.0214 (0.0197), respectively. 


For NRs with SI ranging from 10 to 30 PE (with energy 
approximately from 8 to 20keV nr ), we found that our 
measurement is quite consistent with the results from 
NEST vl.O [23] , thus validating the models used in NEST 
for predicting the response of low energy NRs for a large 
range of drift fields (236 V/cm to 3.93kV/cm) studied in 
this work. 

In this analysis, the excimer-to-ion ratio N ex /Ni for 
ER in Eq. [3] is fixed at 0.06 (which is based on the 
calculation (381 [35]) to be compatible with the treatment 
in NEST [55]. However some measurements indicate 
a larger N ex /N 4 of 0.13T0.07 gD] and 0.20±0.13 gj. 
Additionally in our data, we performed a similar y 2 
analysis as illustrated in Fig. [5j but with both the 4^/N^ 
and N ex /Ni treated as free parameters. The best-fit 
N e;c /Ni obtained for all scanned fields are shown in Fig. [8] 
The mean N ex /Ni is 0.11 with the variance of 0.07. This 
is in consistent with the measurements mm and the 
calculation [381 [351 . 

According to NEST [25], the Thomas-Imel box model 
is not valid for describing the recombination process in 
LXe for ER energy larger than 15keV ee . In this work, 
we also calculated the model-independent photon yields 
at all scanned fields based on the combined energy E c , 


which is defined as: 

SI S9 

Ec = ^PDE + EAF^ Wq ' ^ 

The mean photon yield, P y , is the average number of 
photons generated in LXe per keV energy. Our measured 
photon yields at all scanned fields are shown in Fig. [9] 
along with the photon yields from LUX [87] at 180 V/cm. 
The measured photon yields are observed to deviate from 
the box model when the combined energy is larger than 
8keV ee - Above 8keV ee , the measured photon yields are 
lower than the box model prediction, indicating that the 
Doke-Birk recombination g2] starts to contribute to the 
process. We report the results here and leave the data 
fitting for different recombination models at all energies 
for future publications, together within the NEST group. 

V. CONCLUSION 

In summary, we performed new measurements of 
scintillation and ionization of LXe for low energy 
electronic and NRs at drift fields from 236 V/cm to 
3.93 kV/cm using a three-dimensional sensitive LXe time 
projection chamber. The three-dimensional sensitivity 
allows the removal of edge events which reduces the 
systematic errors. The responses to NRs from our 
measurement are quite consistent with the NEST model. 
But the responses to ERs from our measurement differ 
from the parameters in NEST. 

By using a simulation taking into account the detector 
parameters and all statistical effects, we are able to 
reproduce the electronic and NR bands in Logio(S2/Sl) 
over SI space and obtain the recombination parameters 
for ER at all drift fields using a minimum y 2 method. 
Our obtained recombination parameters for ER bands 
are well fit by the Thomas-Imel box model (SI from 8 to 
40PE), with 4£/Ni values about 30% lower than those 
in the current NEST model. 

We also provide the model-independent ER photon 
yields as a function of the combined energy at all scanned 
field. The photon yields deviate from the box model 
when the energy is larger than ~8keV ee . These data 
are useful in the global analysis and modeling of the ER 
recombination in the low energy region. 

Due to the lack of precise measurement at fixed 
energy below 10 keV ee previously, our new measurements 
provide a set of best available data to be used to predict 
the response of low energy ERs in LXe. In addition to 
that, we measured the response from very low to high 
drift fields, providing useful data to predict the response 
of LXe at different drift fields for future large LXe dark 
matter detectors. 
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FIG. 9: The photon yields P y as a function of combined energy (blue lines) in this work at all scanned fields. The recent result 
from LUX [37] is plotted in the left-upper panel as well. The uncertainties (blue shadow) include those from PDE and EAF. 
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